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Abstract 



We review recent progress in limit laws for the one-dimensional asymmetric simple exclusion 
process (ASEP) on the integer lattice. The limit laws are expressed in terms of a certain Painleve 
II function. Furthermore, we take this opportunity to give a brief survey of the appearance of 
Painleve functions in statistical physics. 
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"It was a pleasant surprise to me that such special functions actually appeared in concrete 
problems of theoretical physics. . ." Mikio Sato [4]. 

1 Introduction 

The appearance of Painleve functions in the 2D Ising model is well-known [37, 64]. Equally well- 
known is that this provided one impetus for M. Sato, T. Miwa and M. Jimbo [48] to develop 
their theory of holonomic quantum fields which connects the theory of isomondromy preserving 
deformations of linear differential equations with the n-point correlation functions of the 2D Ising 
model. 1 

The general consensus in the field of "exactly solvable models" is that correlation functions 
are expressible in terms of Painleve functions only in models that are free fermion models. More 
precisely, one expects that for the appearance of functions of the Painleve type, it is necessary for 
the underlying model or process to be a determinantal process in the sense of Soshnikov [52]. In 
addition to the 2D Ising model, some notable examples where Painleve functions arise in correlation 
functions include the one-dimensional impenetrable Bose gas [21, 28, 33, 34], the Ising chain in a 
tranverse field [41], the distribution functions of random matrix theory [1, 5, 16, 22, 28, 56, 57, 58], 
Hammersley's growth process [7, 8], corner and polynuclear growth models [9, 24, 29, 42, 43] and 
the totally asymmetric simple exclusion process (TASEP) [12, 29, 44]. Universality theorems in 
random matrix theory have extended the appearance of Painleve functions to a wide class of matrix 
ensembles [13, 17, 18, 19, 51]. 2 In recent developments [3, 45, 46, 47] Painleve II appears in the 
long time asymptotics of explicit formulas for the exact height distribution for the KPZ equation 
[32] with narrow wedge initial condition. 

As just noted, one does not expect Painleve functions to arise in correlation functions in models 
that are exactly solvable in the sense of Baxter [11] but are not free fermion models, e.g. 6- vertex 
model, XXZ quantum spin chain, Baxter's 8-vertex model. Having said that, the universality 
conjecture arising in the theory of phase transitions suggests, for instance, that the scaling limit of a 
large class of ferromagnetic 2D Ising models is the same as that of the Onsager 2D Ising model; and 
hence, Painleve functions are conjectured to appear (in the massive scaling limit) in models outside 
of the class of exactly solvable models. This last statement is substantiated by the developments in 
[3, 45, 46, 47]. 

In this paper we review recent progress [59, 60, 61, 62, 63] on the current fluctuations in the 
asymmetric simple exclusion process (ASEP) on the integer lattice Z [35, 36]. ASEP is in the class 

X A complete account of the SMJ theory can be found in the recent monograph by Palmer [39]. 
2 It is also worth noting that due to the close connection of random matrix theory to multivariate statistical analysis, 
these same distribution functions involving Painleve functions are now routinely used in data analysis [30, 31, 40]. 
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of Bethe Ansatz solvable models [23, 25] but only for certain values of the parameters is ASEP a 
determinantal process [29, 44, 49]. That ASEP is Bethe Ansatz solvable comes as no surprise once 
one realizes that the generator of ASEP is a similarity (not unitary!) transformation of the XXZ- 
quantum spin Hamiltonian [2, 50, 65]. Our main results relate the limiting current fluctuations in 
ASEP for certain initial conditions to the TW distributions i*\ and F 2 of random matrix theory 
[58, 59]. Both F\ and F 2 are expressible in terms of the same Hastings-McLeod solution of Painleve 
II [20, 26], see §4.2. 



2 Master Equation and Bethe Ansatz Solution 



Since its introduction in 1970 by F. Spitzer [53], the asymmetric simple exclusion process (ASEP) 
has attracted considerable attention both in the mathematics and physics literature due to the fact 
it is one of the simplest lattice models describing transport far from equilibrium. Recall [35, 36] that 
the ASEP on the integer lattice Z is a continuous time Markov process r\ t where rft(x) = 1 if x € Z 
is occupied at time t, and rj t {x) = if x is vacant at time t. Particles move on Z according to two 
rules: (1) A particle at x waits an exponential time with parameter one, and then chooses y with 
probability p(x,y); (2) If y is vacant at that time it moves to y, while if y is occupied it remains at 
x. The adjective "simple" refers to the fact that the allowed jumps are only one step to the right, 
p(x, x + 1) = p, or one step to the left, p(x, x — 1) = q = 1 — p. The totally asymmetric simple 
exclusion process (TASEP) allows jumps only to the right (p = 1) or only to the left (p = 0). 3 In 
the mapping from the XXZ quantum spin chain, the anisotropy parameter A of the spin chain is 
related to the hopping probabilities p and q by 




the ferromagnetic regime of the XXZ spin chain. 

We begin with a system of N particles and later take the limit N — > oo. A configuration is 
specified by giving the location of the N particles. We denote by Y = {yi, . . . , yjy} with yi < ■ ■ ■ < 
yjv the initial configuration of the process and write X = {x\, . . . , xn} € "L N . When x\ < ■ ■ ■ < xjv 
then X represents a possible configuration of the system at a later time t. We denote by Py (X; t) the 
probability that the system is in configuration X at time t, given that it was initially in configuration 
Y. 

Given X = {x\, ... ,xn} G Z n we set 

X t = { x li ■ ■ ■ , x i-l, x i + 1,^+1, • • • , x n}, X~ = {x 1} . . . ,Xi-i,Xi - l,X i+ i, . . . ,X N }. 
3 It is TASEP that is a determinantal process. 
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The master equation for a function u on Z N x R + is 

j N 

ju{X-t) = (pu(X-;t)+q U (X+;t)-u(X;t)), (1) 
i=i 

and the boundary conditions are, for i = 1, . . . , N — 1, 

u(xi, . . . , Xj, Xj + 1, . . . , xtv; i) 

= pu(xi, . . . ,Xi,Xi, . . . ,x N ;t) + qu(xi, . . . , x» + 1, Xj + 1, . . . ,x N ;t). (2) 
The initial condition is 

u(X;0)=5y(X) when xi < • • • < xat. (3) 

The basic fact is that if u(X;t) satisfies the master equation, the boundary conditions, and the 
initial condition, then Py(X;t) = u(X;t) when x\ < ■ ■ ■ < xn- This is, of course, one of Bethe's 
basic ideas (see, e.g., [10]): incorporate the interaction (in this case the exclusion property) into the 
boundary conditions (2) of a free particle system (1). 

Recall that an inversion in a permutation a is an ordered pair {a(i),a(j)} in which i < j and 
a(i) > a(j). We define [65] 

and then 

A a = Y\{Sa/3 '■ is an inversion in a}. 

We also set 

e (o=pr 1 +?e-i- 

In the next theorem we shall assume p / 0, so the A a are analytic at zero in all the variables. Here 
and later all differentials d£ incorporate the factor (2-Ki)^ 1 . 

Theorem 2.1. We have 

p Y (x; t )= £ / .../ ^ncr w " leE<e(€,)i ^-^ ^ 

aeS N Cr Cr i 

where C r is a circle centered at zero with radius r so small that all the poles of the integrand lie 
outside C r . 

The proof that Py(X;t) satisfies (1) is immediate and the fact it satisfies the boundary condi- 
tions (2) is exactly the same argument as in the XXZ problem [65]. The difficulty lies in showing 
(5) satisfies the initial condition (3). Observe that the term in (5) corresponding to the identity 
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permutation does satisfy the initial condition. Thus the proof will be complete once one demon- 
strates that the remaining n! — 1 other terms sum to zero at t = 0. This is indeed the case (some 
are individually zero and others cancel in pairs) and the result depends crucially upon the choice 
of the contours C r [59]. For the special case of TASEP, p = 1, it follows from (4) and (5) that the 
right-hand side of (5) can be expressed as a JV x JV determinant as first obtained in [49]. 

We note that unlike the usual applications of Bethe Ansatz, it is not the spectral theory of 
the operator that is of interest but rather the transition probability Py(X;t). Thus there are 
no Bethe equations in our approach; and hence, no issues concerning the completeness of the 
Bethe eigenfunctions. Indeed, there is not even an Ansatz in this approach! We remark that this 
result extends with only minor modifications to the solution VP(xi, . . . , xjv; t) of the time-dependent 
Schrodinger equation with XXZ Hamiltonian where the Xj's denote the location of the N "up spins" 
in a sea of "down spins" on Z. 



3 Marginal Distributions and the Large TV Limit 

We henceforth assume q > p so there is a net drift of particles to the left. Here we consider two 
different initial conditions. The first, called step initial condition, starts with particles located at 
Z + = {1,2,...}. The second initial condition is the step Bernoulli initial condition: each site in 
Z + , independently of the others, is initially occupied with probability p, < p < 1; all other sites 
are initially unoccupied. In each of these cases it makes sense to speak of the position of the mth 
particle from the left at time t, x m (t), and its distribution function ¥(x m (t) < x). It is elementary 
to relate F(x m (t) < x) to the distribution of the total current T at position x at time t, 

T(x,t) := number of particles < x at time t; 

namely, 

P(T(x,t) < m) = 1 - ¥(x m+1 (t) < x). 

For this reason we first concentrate on Py(x m (i) < x) and only at the end translate the results into 
statements concerning T. (The subscript Y denotes the initial configuration.) 

Now for finite Y 

P Y (x m (t) = x) = ^2 Py{xi, ■ ■ ■ ,x m -i,x,x m +i, . . . ,x N ;t). 

Xl<---<X m -l<X<X m+ l<---<X N 

Since the contours C r in (5) have r < 1, the sums over x m +i, . . . , xn can be interchanged with 
the integrations in variables £^L\, m + 1 < j < N , and the geometric series summed. To perform 

the sums over x\, . . . ,x m -\ the contours in the £, Xj ( -\ variables, 1 < j < m — 1, must be deformed 



4 



out beyond the unit circle and then the sums can be interchanged with the integrations. This 
deformation beyond the unit circle can be done in such a way as not to encounter any poles of 
the integrand. However, upon deforming these contours back to C r (after the geometric series are 
summed) one does encounter poles; and one finds some remarkable cancellations: only the residues 
from the poles at & = 1 are nonzero. The result is a sum over all subsets of S of {1,...,N} 
with \S C \ < m whose summands involve \S\ -dimensional integrals with contours C r . 4 However, this 
resulting expression for Py(x m (<) = x) is not so useful for taking the N — > oo limit. 

The next step is to expand the contours to Cr, R S> 1. It is then possible to take the N — > oo 
limit in the resulting expressions. The details [59] are involved and they depend crucially upon 
some algebraic identities which we now state. 



3.1 Three identities 



Let 

Identity #1: 



f(i,j) :=p + q^j - &• 



ru/H*),^)) 



\ sgn(cr) 

J^S ~~ 1 ) (^(1)^(2) -I)-" (6t(1)£<t(2) • • • L(N) ~ 1) 



q N(N-i)/2 rL<j(o &) ^ 



n,-(€i-i) 



Identity #2: For N > m + 1, 



E llfVI' II . 

1 jes c 



N - 1 

m 



N 

(i-IH-) 



(7) 



In (7) the sum runs over all subsets S 1 of {1, . . . , N} with cardinality m, and S c denotes the com- 
plement of S in {1, . . . , N}. Here ^ is a slightly modified r-binomial coefficient, r := p/q, 

N _ N 



p-q 

[N}\ := [N] [N — 1] •••[!], 



N 
rn 



[N}\ 



[m]\[N-m]\ 



ra(N-m) 



N 
m 



where 



v 



is the usual r-binomial coefficient. We define 



first proved the simpler identity 



N 



for m < 0. In proving (7) we 



4 This is Theorem 5.1 in [59]. 
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Identity #3: 



e n 



\S\=m <e 
1 1 jes c 



s & ~ & 



N 
m 



We believe that these identities suggest a deeper mathematical structure that is yet to be discovered. 



3.2 Final expression for F(x m (t) < x) for step and step Bernoulli initial condi- 
tions 



We denote by ¥ p the probability measure for ASEP with step Bernoulli initial conditions. For p = 1 
the measure is ASEP with step initial condition. Let 



C m ,k ■= (-l) m g fe (fc-l) T ™(™-l)/2 r -fcm 



k-1 
m — 1 



Observe that c m ^ = when m > k. 
Theorem 3.1 [59, 63]. Assume q > p, then 



V P (x m (t)<x) = 



fe(fc-l)/2 r fc(fc+l)/2 



k>l 



k\ 



Cf 



Cr 



n 



TT ^ TT Si 

ll^-l + Kl-r) 11 1 



(8) 



The contour Cr, a circle of radius i? S> 1 centered at the origin, is chosen so that all (finite) poles 
of the integrand lie inside the contour. 

We remark that for TASEP, p = 0, the above sum reduces to one term; and this term can be 
shown to be equal to a m x m determinant. 

The final simplification results if we use the identity [60] 

1 \ , -.1,, >lfl_^/oTT (()-^ 



det 



f(hj) J l<i,j<k 



= (-i)*(w)' ( *- 1>/2 fi- 



ll 



f(i,j) 



in (8) and recognize the summand, a A;- dimensional integral, as the coefficient of \ k in the Fredholm 
expansion of det(I — XK p ) where K p acts on functions on Cr by 



Jc R 
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where 



*,«,?)=, r;;' f-: } , .r=£. <9) 

P + - f £ - 1 + p(l - t) ? 



Note that when p = 1, the case of step initial condition, the last factor in K p (£,£') equals one. 
Since the coefficient of X k in the expansion of det(7 — XK p ) is equal to 



this fact together with the r-binomial theorem gives the final result for ¥ p (x m (t) < x). 

Theorem 3.2 [59, 63]. Let F p denote the probability measure for ASEP with step Bernoulli initial 
condition with density p and x m (t) denote the position of the mth particle from the left at time t, 
then 

det(I-XK p ) dX 



V P (x m {t) <x)= \ f- - (10) 

where the contour C is a circle centered at the origin enclosing all the singularities at A = t~ j , 
< j < m — 1 and K p is the integral operator whose kernel is given by (9) . 



4 Limit Theorems 
4.1 KPZ Scaling 

The scaling limit that is of most interest is the KPZ scaling limit [32, 54]. In the terminology here 
this scaling limit is 

TYl 

m — > oo, t — > oo with a = — < 1 fixed. 

As we shall see, the limiting distribution will depend upon the relative sizes of a and p 2 . For the 
moment we concentrate on the cases < a < p 2 and a = p 2 with < p < 1. As in any central 
limit theorem, to obtain a nontrivial limit the x in F p (x m (t) < x) must also be scaled (this too is 
part of KPZ scaling). In anticipation of the theorem we set 

x := at + c 2 t 1/3 s 

where the | is the famous KPZ universality exponent [32, 38] and 

a := -1 + 2v^, c 2 := a^/^l - y^) 2 / 3 . 
The two distribution functions that arise in the KPZ scaling limit are defined in the next section. 
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4.2 Distributions F x and F 2 



The distributions F\ and F 2 can be defined by either their Fredholm determinant representations or 
their representations in terms of a Painleve II function. Here we take the latter route. Let q denote 
the solution to the Painleve II equation 



where Ai(x) is the Airy function. That such a solution exists and is unique was proved by Hastings 
and McCleod [26]. 5 Then we have 



The asymptotics of these distributions as x — > oo is straightforward given the large x asymptotics 
of the Airy function; however, the complete asymptotic expansion as x — > — oo has only recently 
been completed [6]. For high-accuracy numerical evaluation of F\ and F 2 , it turns out that it is 
better to start with their Fredholm determinant representations [15]. 

4.3 Limit Laws 

The asymptotic analysis [61, 63] of the Fredholm determinant in the formula for F p (x m (t) < x) in 
(10) required the development of new methods since the operator K p is not of the usual "integrable 
integral operator" form normally appearing in random matrix theory [14, 27, 57]. The main point 
is that the kernel K p has the same Fredholm determinant as sum of two kernels; one has large norm 
but fixed spectrum and its resolvent can be computed exactly, and the other is better behaved [61]. 

We now state the results of this asymptotic analysis. 

Theorem 4.1 [61, 63]. When < p < q, 7 := q — p, 



q" = xq + 2q 3 



satisfying 



q(x) ~ Ai(x), x — > 00 




(11) 



(12) 




(14) 



(13) 



5 A modern account of Painleve transcendents can be found in the monograph by Fokas, et al. [20]. 



8 



Table 1: The mean (pp), variance (<rl), skewness (Sp) and kurtosis (Kp) of Fp, /3 = 1,2. The 
numbers are courtesy of F. Bornemann and M. Prahofer. 



/? 




4 


Sp 




1 


-1.206 533 574 582 


1.607 781034 581 


0.293 464 52408 


0.165 242 9384 


2 


-1.771086 807411 


0.813194 792 8329 


0.224084 203 610 


0.093 448 0876 



This theorem implies a limit law for the current fluctuations. Define 

v = x/t, ai = {l + v) 2 /4, a 2 = 2" 4 / 3 (l - t> 2 ) 2 / 3 . 

Theorem 4.2. When < p < q, 7 := q — p, 

lim P p f T(^,t/ 7 ) - Ql t ^ \ = 1 _ F2( _ s) _ 1 < v < 2p _ 1 

l im P ( n^th)- ai t <\ =1 _ p r_ s) 2 when v = 2 p-l, p< \. (16) 
t^oo y a2t 1 ' 6 J 

For step initial condition with < a < 1 the limit laws are (13) and (15) [61, 62]. When a > p 2 (or 
v > 2p — 1) the fluctuations are of order t 1 / 2 and the limiting distribution is Gaussian, see [63] for 
details. 

For TASEP, p = 0, with step initial condition the limit law (15) was first proved by Johansson 
[29]. For TASEP with step Bernoulli initial condition the limit laws (15) and (16) were conjectured 
by Prahofer and Spohn [44] and proved recently by Ben Arous and Corwin [12]. The fact that these 
limit laws remain essentially identical (the only change is the factor 7 in the time slot) is a very 
strong statement of KPZ Universality. From the integrable systems perspective, these results are, to 
the best of the authors' knowledge, the first limit laws of Bethe ansatz solvable models (outside the 
class of determinantal models) where the correlation functions are expressible in terms of Painleve 
functions. 



5 Conclusions 



Today Painleve functions occur in many areas of theoretical statistical physics. In the case of 
KPZ fluctuations there are now experiments [38, 55] on stochastically growing interfaces where 
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quantities such as the skewness and the kurtosis of Fp (see Table 1), as well as the distribution 
functions themselves, are compared with experiment. In [55] K. Takeuchi and M. Sano conclude 
that their measurements "...have shown without fitting that the fluctuations of the cluster local 
radius asymptotically obey the Tracy- Widom distribution of the GUE random matrices." 
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DMS-0906387 (first author) and DMS-0854934 (second author). 
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